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ABSTRACT 



We have imaged the extremely high- velocity outflowing gas in CO (2—1) and 
(3—2) associated with the shell-like ultracompact Hn region G5. 89— 0.39 at a 
resolution of ~3" (corresponding to ~4000 AU) with the Submillimeter Array. 
The integrated high- velocity (>45 km s _1 ) CO emission reveals at least three 
blueshifted lobes and two redshifted lobes. These lobes belong to two outflows, 
one oriented N-S, the other NW-SE. The NW-SE outflow is likely identical to 
the previously detected Br7 outflow. Furthermore, these outflow lobes all clearly 
show a Hubble-like kinematic structure. For the first time, we estimate the 
temperature of the outflowing g function of velocity with the large velocity 

gradient calculations. Our results reveal a clear increasing trend of temperature 
with gas velocity. The observational features of the extremely high-velocity gas 
associated with G5. 89— 0.39 qualitatively favor the jet-driven bow shock model. 

Subject headings: Hn regions — ISM: individual objects (G5. 89— 0.39) — ISM: 
jets and outflows — stars: formation 



1. Introduction 



Although molecular outflows have b een commonly identified around both low- and high- 
mass young stellar objects (YSOs) (e.g.. iBachiller &: Tafallalll999l : lArce et al.ll2007| ). it is not 
clear how the bulk of the outflowing gas is accelerated. Past observations show that molecular 
outflows can often be divided into two components — the "classical" less-collimated com- 
ponent with low-velocity and the highly collimated component with extremely high-velocity 
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(50— 150 km s _1 , hereafter EHV) flBachiller fc Tafallalll999t iHirano et alJbooa IPalau et al. 
20061 ). Exploring the physical conditions of the EHV molecular gas in details will be helpful 
for clarifying its role in star formation processes. 

The massive star forming region G5. 89— 0.39 (hereafter G5.89) is associated with ener- 
getic CO outflows with velocities up to ~±70 km s _1 from its systemic v elocity, among th e 
highest gas velocities that have been detected in molecular outflows (e.g., [Choi et al.lll9931 ). 



G5.89 also harbors a shell-like ultracompact (U C) H n region powered by a young O-type sta r 



(hereafter Feldt's star) visible in the near-IR (IWood fc Churchwellll 19891 : iFeldt et al.ll2003h . 
Based on the deduced gas temperature structure, ISu et al. J2009h conclude that Feldt's star 



supplies the majority of energy for heating the dense gas enshrouding the UC H 11 region. 



The outflowing gas associated with G5.89 was first identified by lHarvey Forveille 



( 119881 ) in CO (1—0) with the IRAM 30-m telescope. Subsequent studies reported outflows 
with notably different orientations in various tracers. For example, SiO (5—4) observations 
with the Submil limeter Arra\j] (SM A) resolved a bipolar structure along a position angle 



(P.A.) of -30° flSollins et al. 



made with the SMA ( Hunter et al. 



2004) . On the other hand, sub-arcsecond CO (3—2) images 



20081 ) detected a well-collimated N-S outflo w and tenta- 



tively a compact NW-SE outflow associated with the Br7 outflow revealed by iPuga et al. 
( 120061 ). Both CO outflows demonstrate EHV gas, with velocities >50 km s -1 . Furthermore, 
BIMA observations in CO (1-0) and HCO+ (1-0) identified an outflow along the E-W 



direction ( IWatson et al.l 120071 ). These results imply the existence of multiple YSOs toward 
G5.89. Indeed, five sub-mm dust sources (denoted as SMA1, SMA2, S MA-N, SMA-E, an d 
SMA-S) were located in the close vicinity of the UC Hn region G5.89 ( jHunter et al.l 120081 ). 



Here we present SMA observations of G5.89 in CO (2—1) and (3—2) at resolutions 
of ~3". We perform large velocity gradient (LVG) calculations to explore the physical 
conditions of the EHV outflowing gas. We then discuss the nature of the EHV components 
and their connection to the gas accelerating processes. 



2. Observations and Data Reduction 

The observations were carried out with the SMA on 2008 April 17 and May 2 in the 
band of 230 GHz, and on 2006 July 27 and September 10 in the band of 345 GHz. The 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
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phase center was R.A. = 18 /l 00 m 30.32 s (J2000) and decl. = -24°04'00.50" (J2000). At 230 
GHz band, CO (2—1) was observed with seven antennas in the compact configuration, which 
resulted in projected baselines ranging from about 9 m to 120 m (7 to 90 kA). For the CO 
(3—2) observations at 345 GHz band with also seven antennas in the compact configuration, 
the projected baselines ranged from about 9 m to 70 m (10 to 80 kA). The half-power width 
of the SMA primary beam was ~51" at 230 GHz and ~34" at 345 GHz. The total available 
double-sideband bandwidth was 4 GHz. The spectral resolutions of 230 GHz and 345 GHz 
observations were 0.41 M Hz (0.53 km s _1 ) and 0.81 MHz (0.71 km s _1 ), respectively. See 
Ho. Moran. &: Lol (120041 ) for more complete specifications of the SMA. 



For the 345 GHz band observations, the flux calibrator was Uranus and the bandpass 
calibrators were Uranus (July 27 and September 10), Callisto (July 27), and 3C279 (Septem- 
ber 10). The nearby compact radio sources 1626—298 (S ~ 2.2 Jy in July 27 and 1.8 Jy in 
September 10) and 1924—292 ( S ~ 3.4 Jy in Ju ly 27 and 2.9 Jy in September 10) served 
as complex gain calibrators. See ISu et al.l (120091 ) for the calibration information of the 230 
GHz band observations. For both 230 GHz and 345 GHz band observations, the absolute 
flux density scales were estimated to have an uncertainty of ~20% . We calibrated the data 
using the MIR software package adapted for the SMA from the software package developed 
originally for the OVRO M MA JScoville et al.lll993h . We made maps using the MIRIAD 
package (jSault et al.lll995l ). With robust weighting, the synthesized beam size was about 
3.2" x 2.4" at P.A. of 40.8° for the CO (2-1) maps, and about 3.0" x 2.0" at P.A. of 13.5° 
for the CO (3—2) maps. We smoothed our data to 2.0 km s -1 resolution for the analysis 
presented below. The rms noise level in a 2.0 km s -1 velocity bin is ~35 mJy beam" 1 (~0.10 
K) at 230 GHz and ~90 mJy beam" 1 (~0.15 K) at 345 GHz. 



3. CO (2-1) & (3-2) Emissions 

Figure [1] shows channel maps of the blueshifted and redshifted CO (2—1) emissions 
associated with G5.89. For display purposes, we have smoothed the channel maps to a 
velocity resolution of 20.0 km s -1 . The spatial distributions of CO (2—1) and (3—2) agree 
with each other very well, and the velocity extents of the CO (2—1) and (3—2) outflows 
are comparable. Emissions are detected to velocities Vfi ow ~ 80 km s _1 in the redshifted 
lobes and Vfi ow ~ 160 km s _1 in the blueshifted lobes, where Vfi ow = \Vfi ow — Vi sr \ with Vfi ow 
being the apparent outflow velocity and Vi sr (= +9.0 km s -1 ) the systemic velocity of G5.89. 
The EHV line wi ng detected here is significantly broader than that reported by single-dish 



observations (e.g.. lChoi et al.lll993f ). The difference is due to the relatively poor sensitivities 
of single-dish observations. The structure of the CO gas with Vfi ow < 40 km s -1 is extended 
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and complex. Due to the lack of short-spacing data to recover the extended structure filtered 
out by our SMA observations, in Figure [1] the low- velocity components are excluded. 

The left and right panels of Figure |5] show the integrated blueshifted and redshifted 
emissions of CO (2—1) and (3—2), respectively. For both transitions, the emissions are 
integrated over the line wing, with 50 < Vfi ow < 162 km s -1 for the blueshifted gas and 42 
< Vfi ow < 80 km s^ 1 for the redshifted gas. The outflow morphologies seen in CO (2—1) 
and (3—2) are very similar. Although only EHV line wing is integrated, the morphology 
of the EHV gas is complicated — at least three bluedshifted lobes (denoted as B-N lobe, 
B-NW lobe, and B-S lobe) and two redshifted lobes (denoted as R-N lobe and R-S lobe) 
are revealed in the vicinity of the UC Hn region G5.89. There appears to be an additional 
redshifted lobe located ~4" southeast to Feldt's star. 

Both B-S lobe and R-N lobe have been identified by previous SMA CO (3—2) observa- 
tions at sub-arcsecond resolution, and were proposed to be a p air of the N-S bipolar outflow 
driven by the sub-mm dust source SMA1 Jrlunter et al.ll2008h . The authors also proposed 



a second CO ou tflow along P.A. about —55° associated with the Br7 outflow reported by 



Puga et al.l (120061 ). Our observations confirm the existence of this NW- SE outflow, co n sistin g 



of the B-NW lobe and R-S lobe. The B-N lobe was also identified bv iHunter et all (120081 ). 
although its driving source is less clear. The existence of the B-N lobe could be explained 
by the N-S outflow being close to the plane of the sky, but this is very unlikely from the 
kinematic viewpoint. After correcting for projection, the N-S outflow, if in the plane of sky, 
will have molecular gas with velocity up to >500 km s _1 , about a factor of 3 higher than the 
highest gas velocities that ever identified in molecular outflows. Figure [3] shows the position- 
velocity (P-V) diagrams of the CO (2—1) and (3—2) along P.A. 175° (left panel) centered 
on the dust source SMAl and 130° (right panel) centered on the dust source SMA2. The 
Hubble-like kinematic feature, i.e., the outflow velocity increasing with the distance from 
the central YSO, is clearly discerned in both outflows with a velocity gradient about 1000 
km s^ 1 pc _1 . 



Assuming an inclination of 45° and a distance of 1.28 kpc (IMotogi et al.ll201ll ). the 
inferred dynamical timescales range from 400 years to 650 years. Since the axes of these 
outflows are unlikely close to either the plane of the sky or the line of sight, the derived 
dynamical timescales should be accurate to within a factor of 3—4. The observed outflow 
lobes therefore represent the most recent outflow activities in this region. 
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4. Line Ratios & Excitation Conditions of the Outflowing Gas 

Given the fairly good signal-to-noise ratio of the CO emission detected toward G5.89 
even at EHV line wing, it is feasible to estimate the physical parameters with a relative nar- 
row velocity bin (10 km s" 1 ) and investigate their velocity-dependence. Table [1] summarizes 
the observed CO (2—1) intensities (in brightness temperature scale) as well as the intensity 
ratios of CO (3—2) and (2—1) as a function of gas velocity for all the observed lobes. 

For a proper comparison, both CO (2—1) and (3—2) maps were reconstructed with a 
clean beam of 3.4" in order to match the resolutions. For each velocity bin, furthermore, we 
simply report the value calculated from the CO (2—1) peak toward each lobe rather than 
the value averaged over entire lobe to avoid the difficulty of emission separation from various 
lobes, in particular, at relatively low-velocity channels. We emphasize that for each outflow 
lobe, the peak positions of CO (2—1) and (3—2) in various velocity bins agree very well with 
each other, with both median and mean offsets of ~0.25", and no systematic offsets can be 
discerned. Consequently we do not expect any systematic bias for the physical conditions 
estimated afterward. As listed in Table [IJ most deduced line ratios are in the range of 
0.6-1.2. 

To estimate the physical conditions of the outflowing gas, we performed the LVG analysis 
with the code wri tten by L. G. Mundy and implemented as part of the MIRIAD package 



( ISault et al.l Il995l ). We assumed the canonical CO fractional abundance of 10~ 4 and a 
velocity gradient of 1000 km s" 1 pc _1 as estimated from the Hubble-like kinematic structure 
(in Figure [3]). Figure H] shows the gas temperature estimated from the LVG calculations 
versus gas velocity. The plot exhibits a clear increasing trend of gas temperature with 
outflow velocity. Such a trend can be discerned toward all the lobes, although the trend of 
the B-N lobe is less clear. In the case of B-NW lobe, there appears to be a temperature 
jump at outflow velocities of ~150 km s _1 , with gas temperature ranging from 1200—1800 
K for Vfi^ >150 km s _1 and <150 K for Vfi ow <150 km s -1 . The estimated temperature 
of the above highest-velocity bins in the B-NW lobe will be reduced to 400—500 K hence 
eliminating the above-mentioned temperature jump if the assumed CO fractional abundance 
is lowered by a factor of 2. The occurrence of abundance reduction may be attributed to 
partial dissociation or even ionization of CO molecules under high-excitation conditions. 

A few factors may contribute to the uncertainty of the temperature estimates reported 
above. Given the flux calibration uncertainty of 20%, the gas temperature obtained with LVG 
calculations would vary by a factor of 10 or so. In Figure 4, we show with the blue-shaded 
area the 1-a temperature range when allowing flux variations up to 20% in the B-NW lobe. 
Although the temperature uncertainty appears fairly large, we emphasize that the general 
increasing trend of gas temperature with outflow velocity is indeed robust. For a given flux 
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calibration error, the line ratios will be affected in the same way in all channels and render 
a simple temperature shift across all velocity bins. For the LVG calculations, variations in 
input parameters also lead to variations in the derived temperature by a similar factor and 
hence do not change the general increasing trend on temperature with gas velocity. For 
example, the derived temperature will decrease by about a factor of 3 if the adopted CO 
fractional abundance is reduced by half or the velocity gradient increases twice. Although 
variation in velocity gradient or CO abundance across the line wings could mimic the derived 
temperature trend with velocity, the Hubble-like kinematic feature shown in Figure [3] makes 
this confusion quite unlikely. A precise gas temperature estimate will need observations with 
better flux calibration. Alternatively, observations in higher transitions of CO will also help 
to constrain hot gas temperature. 



5. The Origin of the EHV Components and Their Connection to the Gas 

Accelerating Processes 

Observations have shown that the kinetic temperature of high-velocity molecular gas in 
outflows associated with YSOs can be as high as 100—1000 K. For example, the i nner SiO 
knots in HH211 jet have a temperature excess of 300—500 K (IHirano et al.l 120061 ). In the 
case of extremely active molecular jets in L1448-mm, the kinetic tempe rature of the EHV 
bullets close to the YSO is estimated to be > 500 K (INisini et al.l 120071 ) . Gas temperature 



estimated from near-IR Ho observations typically ranges from 1000—2500 K ([Richer et al. 



2000l ; iDionatos et al.l 120101 ) . Both hot H 2 emission and warm SiO bullets are thought to 
trace shocked molecular gas. The broad velocity dispersions revealed from the kinematic 



2001; 


Hirano et al. 


2006; 


Su et al. 


2007) 



Could the heating of the EHV gas associated with G5.89 be dominated by shocks? 
Indeed, shock indi cators have been detected around the highest velocity gas of the N-S and 



NW- SE outflows flHofner fc Churchwelllll996l : iKurtz. Hofner. fc Alvarezll2004l : iHunter et al. 



2008 ; Puga et al. 20061 ) . For example, as shown in Figure [2] the highest- velocity gas in the 
R-N lobe is assoc iated with masers of H2O and Class I CH3OH as well as H2 knots (labeled 
as CI and C2 by iPuga et al.l (120061 )). Furthermore, the positions of the highest velocity gas 
detected in the B-NW and B-S lobes coincide with the H 2 knots B and group A, respectively, 
very well. The highest velocity in the B-NW lobe is also associated with masers of NH 3 and 
Class I CH3OH. The hot (about 150-1800 K) molecular gas in the highest-velocity bins of the 
N-S and NW-SE outflows can be naturally interpreted as shock-heated gas, while the lack 
of velocity spread in the outflow tips (Figure [3]) can be attributed to insufficient sensitivity. 
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It is expected that the outflow driving process provides energetics to not only acceler- 
ate but also heat the (ambient) gas. Estimations made from a simple jet-driven bow shock 
model predict temperature risi ng with outflow velocity and distance from the driving source 
f lHatchell. Fuller fe Laddlll999h . whi le such trends are different f r om th e predictions of other 
classes of models (lArce et al.l 120071 ). lHatchell. Fuller fe Laddl (119991 ) further argued that 
comparing with the molecular cooling mechanisms, the heating (as a consequence of accel- 
eration) is sufficient to maintain outflowing gas temperature a few times higher than that of 
ambient materials. Since the Hubble-like kinematic structure can also be reproduced by the 
jet-driven bow shock model (ILee et al.ll200ll ; lArce et al.ll2007l ). all above-mentioned features 
of the EHV gas associated with G5.89 can be qualitatively interpreted by the jet-driven 
bow shock model. We note that most, if not all, available bow shock models have parame- 
ters typical of low-mass outflows. Models with physical conditions more similar to outflows 
from high-mass young stars will be desired to have quantitative comparisons between ob- 
servational results and model predictions. Observations with better sensitivities and spatial 
resolutions are also essential to search for the bow shock signatures in both morphology and 
kinematics. 
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Fig. 1. — Channel maps of G5. 89— 0.39 in CO (2—1) emission at high velocity. The local 
standard of rest velocity of each channel's center is indicated in the upper right corner of 
each panel. The LSR velocity of the system is about +9.0 km s -1 . Solid contours are at 4, 
10, 20, 40, 60, 100, 150, 250, 350, and 450 x a, and dotted contours indicate -4, -10, -20, 
—40, and —60 x a, where a — 11 mJy beam -1 . The dotted circle shown in the first panel 
represents the primary beam of the SMA observations. The synthesized beam is shown in 
the lower left corner of the first panel. The asteri sk represents the position of the H n region 
powering star, i.e., Feldt's star (IFeldt et al.l 120031 ). 
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Fig. 2. — Extremely high-velocity molecular outflows of G5. 89— 0.39 in CO (2—1) (left) 
and (3—2) (right) imaged with the SMA. For both transitions, the blueshifted and red- 
shifted lobes are integrated over the line wing, with 50 < Vfi ow < 162 km s" 1 for the 
blueshifted gas and 42 < Vfi ow < 80 km s -1 for the redshifted gas. In the left panel, 
the overlaid color scales represent the f ree-free emission a t 2 cm (ITang et al.ll2009i ). filled 
trian gles are the n e ar-IR H2 knots from iPuga et al.l (120061 ). filled square is the NH3 maser 
from lHunter et ajj f|20081), open circles are Class I methanol maser positions (components 1 



and 2 from 



Kurtz. Hofner. fc Alvarez! ( 2004 )) and open s quares are positions of water masers 



the crosses mark 



(components 1 and 3 from Hofner fc Churchwell ( 19961 )). In each pane 
the position of the sub-mm dust co mponents reporte d by iHunter et al.l (120081 ). the asterisk 



3DC 

marks the posi tion of Feldt's st ar (iFeldt et al.l 120031 ) . and the filled circle labels the Br7 
outflow center (IPuga et al.l 120061 ) . The synthesized beam is shown in the lower left corner. 
The two black lines represent the axes of the position- velocity diagrams plotted in Figure [3J 
In the left panel, contours are 4, 10, 20, 40, 80, 120, 160, 200, and 240 x a, where a = 0.53 
and 0.31 Jy beam -1 km s _1 , respectively, for the blue- and redshifted components. In the 
right panel, contours are 4, 10, 20, 40, 80, 120, 160, and 200 x cr, where a = 1.36 and 0.80 
Jy beam -1 km s _1 , respectively, for the blueshifted and redshifted components. 



- 12 - 




Fig. 3. — Position- velocity diagram of G5. 89— 0.39 in CO (2—1) (contours) and (3—2) 
(grayscales) along the P.A. of 175° centered at SMAl (left) and P.A. of 130° centered at 
SMA2 (right). In each panel, the horizontal line marks the LSR velocity of the system (i.e., 
+9.0 km s _1 ), and the vertical line indicates the position of SMAl (left) and SMA2 (right). 
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Fig. 4. — Temperature of the outflowing gas deduced from the LVG calculations versus the 
corresponding gas velocity. An increasing trend of gas temperature with outflow velocity 
can be discerned toward all lobes. The blue-shaded area represents the 1-cr temperature 
uncertainty when allowing flux variations up to 20% in the B-NW lobe. Note that the 
outflow velocity Vf\ ow shown here is related to the velocity Vfi ow listed in Table 1 by the 
relation: Vfi ow = \Vfi ow — Vi sr \ with Vi sr = +9.0 km s" 1 the systemic velocity of G5. 89— 0.39. 
The lobes of N-S and NW-SE outflows are represented by triangles and circles, respectively. 
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-100 to -90 


1.78 


0.73 


1.78 


1.10 
















-90 to -80 


2.04 


0.78 


2.29 


1.07 


0.81 


0.60 












-80 to -70 


2.84 


0.82 


1.92 


1.15 


3.70 


0.92 












-70 to -60 


3.99 


0.94 


2.14 


1.13 


6.61 


0.91 












-60 to -50 


6.35 


0.94 


2.61 


1.19 


10.54 


0.87 












-50 to -40 


11.08 d 


0.87 d 


11.08 d 


0.87 d 


14.41 


0.84 












-40 to -30 


18.92 6 


0.84 e 


18.92 e 


0.84 e 


18.92 6 


0.84 e 













a measured at the peak position of the CO (2—1) map 

b 3— 2/2— 1 measured at the peak position of the CO (2—1) map 

c CO (3—2) blended with 34 S02, which can also be discerned from Figure [3] 

including B-NW lobe and B-N lobe 

"including B-NW, B-N, and B-S lobes 

f including R-N lobe and R-S lobe 



